Non-universal bound states of two identical heavy fermions and one light particle 
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We study the behavior of the bound state energy of a system consisting of two identical heavy 
fermions of mass M and a light particle of mass m. The heavy fermions interact with the light particle 
through a short-range two-body potential with positive s-wave scattering length a s . We impose a 
short-range boundary condition on the logarithmic derivative of the hyperradial wavefunction and 
show that, in the regime where Efimov states are absent, a non-universal three-body state "cuts 
through" the universal three-body states previously described by Kartavtsev and Malykh [O. I. 
Kartavtsev and A. V. Malykh, J. Phys. B 40, 1429 (2007)]. The presence of the non-universal 
state alters the behavior of the universal states in certain regions of the parameter space. We show 
that the existence of the non-universal state is predicted accurately by a simple quantum defect 
theory model that utilizes hyperspherical coordinates. An empirical two-state model is employed to 
quantify the coupling of the non-universal state to the universal states. 



I. INTRODUCTION 

The unprecedented control over ultracold atomic Fermi 
gases in optical lattices has made them prime candidates 
for studying and engineering novel quantum phases as 
well as probing fundamental theories such as the BEC- 
BCS crossover [THB]. This progress has been facilitated 
by the tunability of the two-body interactions using Fes- 
hbach resonances or by changing the lattice confinement. 
If the interspecies two-body s-wave scattering length a s is 
much larger than the range r of the underlying two-body 
potential, the few- and many-body behavior of equal- 
mass two-component Fermi gases is universal, i.e. com- 
pletely determined by a s . 

Presently significant efforts are directed at creating 
ultracold atomic Fermi gas mixtures composed of two 
chemically distinct species [THS]- This introduces a new 
parameter, the mass ratio k between the two species. 
This new parameter affects the many-body physics of 
the system, allowing one to realize novel quantum phases 
such as the interior gap superfiuid jTUj- Here we show 
that at the few-body level this additional degree of free- 
dom leads to new three-body resonances which may 
destabilize the system, making it harder for experiments 
to explore novel quantum phases with unequal-mass mix- 
tures. In particular, we study a system of two identical 
heavy fermions with mass M, which interact with a light 
particle through a short-range potential with positive s- 
wave scattering length a s . 

Previous studies revealed two intriguing properties. 
First, Kartavtsev and Malykh [TT] predicted the exis- 
tence of a universal trimer state with energy _E u .i for 
Ki < k < K2 and the existence of two universal states 
with energies E u i and E u 2 for «2 < k < 13.606; K\ 
and K2 were found to be 8.173 and 12.917, respectively. 
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Second, Endo et al. [12] investigated how the universal 
trimer states, which are completely determined by the 
s-wave scattering length and the mass ratio k, are con- 
nected to Efimov trimers, which have been predicted to 
exist for k > 13.606. By analyzing the trimer system 
within the framework of the Skorniakov-Ter-Martirosian 
equation with a momentum cutoff A c , Endo et al. pre- 
dicted the existence of a third class of trimer states, 
termed crossover trimers, which were shown to contin- 
uously connect the universal trimers described by Kar- 
tavtsev and Malykh and Efimov trimers. 

Our study employs, as in Ref. [TT], the hyperspherical 
coordinates. However, while Ref. [TT] enforced that the 
hyperradial wavefunction vanishes at hyperradius R = 0, 
we explore the entirety of physically allowed boundary 
conditions by introducing a short-range three-body or hy- 
perradial phase S(Rq). We determine the eigenspectrum 
as a function of the value of the three-body phase 8, the 
hyperradius Rq at which the hyperradial boundary condi- 
tion is imposed and the mass ratio K. The universal states 
of Kartavtsev and Malykh are recovered for Rq — > and 
S(Rq) = 7r/2. However, for other boundary conditions 
we find deviations from universality, which are linked to 
the existence of a non-universal three-body state. Anal- 
ogous non-universal three-body states have previously 
been shown to exist [Jj2 [2] (see also Refs. [T3HT7] ) in the 
a s — > oo limit. The existence of the non-universal state 
for positive a s is described accurately within a quantum 
defect theory (QDT) framework. Moreover, within a two- 
state model, deviations from universality are explained 
as being due to the coupling between the non-universal 
state and the universal states. Our work provides a sim- 
ple intuitive Schrodinger equation based description of 
the energy spectra of heavy-light trimers and an alter- 
native means to understanding the connection between 
universal trimers and Efimov trimers. 

The remainder of this paper is organized as follows: 
Sect ion [H] describes the hyperspherical framework. Sec- 
tion |III| determines the three-body energies as functions 
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of S(Ro) and k by solving the hyperradial Schrodinger 
equation numerically. Section |IV| develops an analytical 
description, which accounts for the universal and non- 
universal states of the energy spectrum. Finally, Sec. |V| 
concludes. 



II. SYSTEM HAMILTONIAN 

We model the interactions between the heavy and light 
particles by a zero-range two-body pseudopotential with 
s-wave scattering length a s . If we denote the heavy 
fermions of mass M as 1 and 2 and the light particle 
of mass to as 3, the Hamiltonian is given by 



Hfnt = 



2M 



(V 2 +V 2 ) 



2 TO 
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where Fj is the position vector of the ith particle, V 2 is 
the Laplacian of the ith particle and Vint = V ps (fi3) + 
V ps (f 2 3) with 



V pB (nj) = 



2nH 2 a s 

H2b 



d 
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Here \x 2 \, is the reduced mass of the heavy-light pair, 
A^2b = ™^pr, an d r VJ = \ fij\ = |f, - fj\. The pseudopo- 
tential V ps imposes the Bcthe-Peierls boundary condition 
on the three-body wave function in the limit r,-j — > 0. 

To solve the Schrodinger equation for the Hamilto- 
nian -fftot, we separate off the center of mass motion 
and write the relative wave function \& in terms of the 
hyperradius R and five hyperangles, collectively denoted 
by SI [TTJ [18]. The hyperrradius R, which provides a 
measure of the overall size of the system, is defined by 
^R 2 = M[(n - R cm ) 2 + (f 2 - R cm ) 2 } + m(r 3 - i? cm ) 2 , 
where fi is the three-body reduced mass associated with 
the hyperradius, fi = to,k/\/2k + 1, and R cm denotes the 
center of mass vector. We expand the relative wave- 
function in terms of a set of weight functions F n (R) 



and adiabatic channel functions <f> n (R; SI), which depend 
parametrically on the hyperradius R [19] . 



vK(i?, fi) = J2 R- 5/2 F n (R)<£ n (R; SI). 



(3) 



The adiabatic channel functions SI) satisfy the hy- 

perangular Schrodinger equation at fixed R, 



A 



2^R 2 



V int (R, SI) $„(i?; SI) = [s 2 (R) - 4] <S> n (R; SI). 

In Eq. Q, A denotes the grand angular momentum 
operator, which accounts for the kinetic energy associ- 
ated with the hyperangles SI [15]. Inserting Eq. ^ into 
the relative Schrodinger equation yields a set of coupled 
equations for the weight functions F n (R). 

In the following, we employ the adiabatic approxima- 
tion, which neglects the coupling between the different 
channels 18, 20, 21!. This approximation has been shown 
to provide a qualitatively correct description for a num- 
ber of three-body systems [TTJ |2"2"H2"4"] . In this approxima- 
tion, the hyperradial Schrodinger equation for the lowest 
adiabatic channel reads 
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F (R) = EF (R), 
(5) 

where Qoo(R) is the diagonal correction to the adiabatic 
energies, 

Qo (R) = (^o(R\Sl)\-^ 2 W^))- (6) 

We determine the three-body energies E using a two 
step process. First we find the hyperangular eigenval- 
ues sq(R) and coupling elements Qoo(R). Then we solve 
the radial Schrodinger equation, Eq. ([5]). For states 
with L 7 ' = 1~ symmetry and zero-range interactions, 
the scaled hyperangular eigenvalues sq(R) can be ob- 
tained semi-analytically by solving the transcendental 
equation [Til IT5] 



R 



{s 2 - 1)(1 + 2 K ) 1/4 sec(^f ) [ SoK 2Fl (i( 3 _ So ) ; 1(3 + So ), 



3(1 + k) sm(^f) 



3s (l+Ac)3/2 
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The hyperangular eigenvalues are completely determined 
by R/a s , k and L n . In the limit R/a s — > 0, Sq goes to a 
constant. For the purpose of this study we are only con- 
cerned with positive values of Sq(R) in the small R/a s 
limit, i.e. we only consider mass ratios k < 13.606. Fig- 
ure [l] shows the scaled eigenvalue Sq(Rq) as a function 
of k for Ro/a s — 0.0001. As we discuss in more detail 
below, the value of sq(R) determines the allowed short- 
range boundary condition of F (R). In the large R/a s 
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limit, Qqq vanishes and the quantity - ^° 2 ^2 ap- 
proaches —E 2 b, where E 2b denotes the dimer binding 
energy, E 2h = h 2 /(2fi 2h a 2 s ). 

Equation ^ has two linearly independent solutions, 
which scale as R s o+ 1 / 2 anc l ^-so+1/2 m tne sma n R/a s 

limit and are referred to as the regular and irregular so- 
lutions, respectively. For sq(R) > 1 and Ro/a s — > 0, 
the irregular solution is not normalizable and thus does 
not contribute, implying that Fq(R) goes to zero as 
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FIG. 1. The hyperangular eigenvalue so(Ro) as a function 
of the mass ratio k for Ro/a s — 0.0001. The dashed lines 
labeled I and II mark the mass ratios Ki and K2, at which 
the universal states with energies E u ^ and first become 
bound. 



R — » [15HT7] . On the other hand, in the regime 
< so < 1 both solutions are well-behaved and must 
be included when constructing the general solution [T5T - 
117] . We parameterize the short-range boundary condi- 
tion of the wavefunction F Q (R) at R using the logarith- 
mic derivative £(F (R )), 
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FIG. 2. (Color online) The dotted, dashed, and dash-dotted 
lines show the scaled three-body energy E/E2b as a function 
of 5 using Ro/a s = 0.0001, 0.0003 and 0.001, respectively, for 
(a) k = 8.6 and (b) k — 10. The solid horizontal lines show the 
energy E n> x of the universal state and the solid vertical lines 
show the energy E nu of the non-universal state (see Sec. IV I 
for R /a s = 0.001. 



where — 7r/2 < 6(Rq) < ir/2. Using this parameteriza- 
tion, we cover all possible short-range phases. In the 
special case of S(Rq) = ir/2 and R/a s — > 0, this bound- 
ary condition and the resulting three-body energies agree 
with those of Ref. [TT] . We refer to this boundary condi- 
tion as "universal boundary condition" since the three- 
body states are completely determined by the regular so- 
lution. We refer to the corresponding states as "universal 
states" with energies E n ^\ for n\ < k < K2, and E u ^i and 
-^11,2 for K2 < K < 13.606. In the next section we discuss 
deviations from universality that occur if 8(Ro) is not 
equal to 7r/2. These deviations increase with increasing 
k (for a fixed i?o) and increasing Rq (for a fixed k). 



III. NUMERICAL TREATMENT 

In this section we examine the behavior of the three- 
body bound states with L u = 1~ symmetry, obtained 
numerically using the shooting algorithm, as a function 
of the logarithmic derivative boundary condition for se- 
lected mass ratios. The three-body energies E for k = 8.6 
and 10 are shown in Figs. |2ja) and (b), respectively. 
Dotted, dashed and dash-dotted lines show E as & func- 
tion of the three-body phase d for i?o/a s = 0.0001, 0.0003 
and 0.001, respectively. In Fig.[2]the three-body energies 



have been scaled by the zero-range two-body binding en- 
ergy i?2b- The three-body state becomes unbound with 
respect to the break-up into a dimer and an atom for 
E/E 2 b > — 1. For k = 8.6 and 6 = 7r/2, the system sup- 
ports one three-body bound state. For Ro/a s = 0.0001 
and k = 8.6, E/E^b is nearly constant for 6 > —0.5. At 
5 m —0.5, referred to as the critical angle 8 c (Ro), the 
energy rapidly goes to a large negative value and a sec- 
ond bound state, whose energy is approximately equal 
to £u,i, is supported for 6 < 6 c (Rq). We refer to the 
feature in the vicinity of 8 C as "three-body resonance". 
As Ro/a s increases [see dashed and dash-dotted lines 
in Fig. [2]ja)] the width of the three-body resonance in- 
creases. Note, however, that the deviations from E u ,i 
are small for all Ro/a s considered, except for three-body 
phases very close to 5 c (Rq). As k increases [see Fig.[2jb)], 
the overall behavior of the energy spectrum is unchanged. 
The key trends with increasing k are that, at a fixed Rq, 
the energy away from <5 c (i? ) becomes more negative and 
both 6 c (Rq) and the "width" of the three-body resonance 
increase (see also symbols in Figs. [5] and [7]). 

For k> K 2 , the three-body system supports a second 
bound state. As an example, Fig.|3]shows the three-body 
bound state energy for n = 13 and Ro/a s = 0.0001 as a 
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function of S. Away from <5 c (i?o) ~ 0.15, there exist two 
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FIG. 3. The solid lines show the scaled three-body energies 
E/E 2 b as a function of S for k — 13 and Ro/a s — 0.0001. The 
three-body resonance is located at S c m 0.15. 
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4. The solid line shows the wave function Fq(R) for 
k = 10, R /a s = 0.0001, E/E 2h = -1.4, and 6 = -0.33. 
As can be seen in Fig. [2jb) , this angle is close to S c . For 
comparison, the inset shows the universal wavefunction for 
S — 1.53 (i.e. away from the three-body resonance). 



bound states whose energies depend weakly on S. For 
S > S c , the corresponding hyperradial functions Fq(R) 
possess nodes and 1 node, respectively. For S < 6 C , 
the corresponding hyperradial functions Fq(R) possess 1 
node and 2 nodes, respectively. This reflects the fact that 
a new bound state is being pulled in at 6 ps S c (Rq). 

In general, we find that S c (Ro) depends relatively 
weakly on i?o as long as Ro/a s <C 1. Moreover, the 
three-body energy depends relatively strongly on Rq in 
the vicinity of S c but comparatively weakly on R away 
from S c . This suggests that the states near S c and away 
from 5 C can be classified as non-universal and univer- 
sal, respectively. This interpretation is corroborated by 
our analysis of the hyperradial wavefunction F (R) for 
k = 10, Ro/a s = 0.0001 and S w 5 C . The main part 
of Fig. [4] shows that the wavefunction has an apprecia- 
ble amplitude in the small R/a s region, signaling non- 
universal behavior. Away from S c , in contrast, Fq(R) 
has a vanishingly small amplitude in the small R/a s re- 
gion (see inset of Fig. [4]). We find that F (R) depends 
fairly weakly on R when 6 is away from 6 C , lending fur- 
ther support to our assertion that the three-body system 
behaves universally in this region. For larger Ro/a Sl we 
observe deviations from universality for a larger range of 
three-body phases. Similar non-universal behavior has 
been previously reported for the three-body system at 
unitarity [131 114] and for the four-body system with pos- 
itive s-wave scattering length |25j . 



IV. ANALYTICAL TREATMENT 

Section |IV|\ applies a QDT framework to predict the 
short-range phase S c (Rq) at which the non-universal state 
first becomes bound. Section |IVp develops a two-state 
model to describe the behavior of the universal and non- 
universal states as a function of 5(Rq). 



A. QDT Treatment 

To predict the short-range phase S c (Rq) at which the 
non-universal state first becomes bound with respect 
to the break-up into a dimer and an atom we apply 
QDT [HI [TJ 023 [27]. In the short-range limit, the wave- 
function Fq(R) can be approximated by [13l fT4] 

Fsn(R) = 

^R[Js {R){^§-R) ~ tan(Tra) Y So{R) {sJ^R)l{S) 

where J So (R) and Y So ^ denote the Bessel functions of 
the first and second kind, respectively. The quantum de- 
fect a controls the relative contribution of the regular 
solution J So (r) and the irregular solution Y So ^. A new 
three-body state is expected to be pulled in when the hy- 
perradial solution is dominated by the irregular solution, 
i.e. for a = 1/2. The critical angle S c (Rq) is then given 
by 

tan«5 c (i? ) = Re[i? £(F SR (i? ))], (10) 

where Re denotes the real part and Fsn(Ro) is evaluated 
for E = ~E 2 b and a = 1/2. 

The solid, dotted and dashed lines in Fig. [5] show 
8 c (Rq) determined using the QDT framework for 
R a /a s = 0.0001, 0.0003 and 0.001, respectively. It is 
interesting to note that <5 c (i?o) is very weakly dependent 
on Rq for small k. The symbols in Fig.[5]show the critical 
angle S c (Ro) obtained by analyzing the numerical solu- 
tions of the hyperradial Schrodinger equation, Eq. ([H]), 
for Ro/a s = 0.0001. The numerical results are in excel- 
lent agreement with the QDT prediction. For k > 13 
(see inset of Fig. [5]) , the dependence of 5 C on i? becomes 
more pronounced. Figure [5] shows that the three-body 
system supports a non-universal bound state not only 




FIG. 5. (Color online) The critical short-range phase 5 c (Ro) 
as a function of k. The solid, dotted and dashed lines show 
the QDT prediction for R /a s = 0.0001, 0.0003 and 0.001, re- 
spectively. The circles show <5 c (i?o) obtained by analyzing the 
numerical solutions of the hyperradial Schrodinger equation, 
Eq. |5|, for Ro/a s = 0.0001. The inset shows a blow up of 
the large n region. 



-0.01 



-0.02 



1.5 -1.0 



-0.5 0.0 0.5 

SSr 



1.0 1.5 



FIG. 6. (Color online) Thick solid, dotted and dashed lines 
show the quantity (E — £ , u ,i)/i?2b as a function of 8 — S c for 
Ro/a s = 0.0001 and k = 10, 12.4, and 13, respectively. The 
circles show the quantity (E — E^^)/ ' Eih for Ro/a s — 0.0001 
and k = 13. The thin dashed lines show the results of the 
two-state model. 



for k > 8.619, but also for k < 8.619 [for k = 8.619 
one has sq(0) = 1]. Thus, it may be surprising at first 
sight that the system supports a non-universal state for 
k < 8.619 since the irregular solution cannot be normal- 
ized if so > 1 f° r Ro/o- s — > 0. However, since we im- 
pose the boundary condition at a finite i?o /ci s and not at 
Ro/ a s — 0, the resulting wavefunction can, even though 
it contains an admixture of the irregular solution, be nor- 
malized. Correspondingly, non-universal states can exist 
for k < 8.619. In fact, we find that the system supports a 
non-universal bound state even for k < k±, i.e. for mass 
ratios where universal states are not supported. 



B. Two-state model 

In this section we develop a two-state model that de- 
scribes the behavior of the universal and non-universal 
states as a function of 6(Rq). In our model the three- 
body resonance is an avoided crossing at 5 c (Rq) between 
the universal state with energy E u and a non-universal 
state with energy £" nu . In Fig. [6] we plot the quantity 
(E — E u )/E2b as a function of S — S c . The thick solid, 
dotted and dashed lines correspond to k = 10, 12.4 and 
13, respectively, with Ro/a s = 0.0001 and E u = I? u ,i- 
The circles correspond to k = 13, Ro/a s — 0.0001 and 
E u = J5 U ,2- 

Motivated by Fig. |6j we write the energy of the non- 
universal state as E nu = E u + A 9 h R -2 (5 — <5 c (-Ro)), where 
A is a dimensionless scaling constant; in our analysis we 
use A = 0.00007. The term proportional to A determines 
the slope at which the non-universal state crosses the 
universal state. In our parameterization of E nu , 5 c (Rq) 
should be interpreted as the three-body phase at which 
the universal and non-universal states cross. Our numer- 



ical results show that this crossing point is nearly iden- 
tical to the QDT three-body phase at which the non- 
universal state first becomes bound. Hence we use the 
QDT value in our two-state model. The dimensionless 
two-state Hamiltonian is then 



H 



Eu/E 2 b 



P/E. 
E nu /E. 



2b 
2b 



(11) 



where /3 = /3(Rq,k) is the coupling between the univer- 
sal and non-universal states. Within this model, /3 is 
real and the scaled three-body energies E/E 2 b are given 
by the eigenvalues of the two-state Hamiltonian. We ap- 
ply the two-state model separately to the universal states 
with energies E Uy i and E Uy2 , yielding separate (3 values 
for the two universal states. The values of (3 are ob- 
tained by fitting the energy spectrum predicted by the 
two-state Hamiltonian to the numerically determined en- 
ergies. Thin dashed lines in Fig.|6]show the results for the 
two-state model. It can be seen that the two-state model 
provides a quantitatively correct description of the three- 
body spectra. Moreover, it provides an intuitive physical 
picture in which deviations from universality arise due to 
the coupling of the universal states to the non-universal 
state. 

The circles, diamonds and squares in Fig. [7] show j3 for 
E n = E uA and R /a s = 0.0001,0.0003 and 0.001, respec- 
tively, as a function of k. The value of (3/ E 2 ^ increases 
with increasing n for fixed Rq. Moreover, the dependence 
of (3 / 1 E 2 b on the value of Rq increases with increasing 
k. Since (3 determines the coupling between the univer- 
sal state and the non-universal state, it can be used to 
quantify the deviations from universality. For mass ra- 
tios larger than those considered in this work, the heavy- 
light trimer system supports three-body bound states 
with Efimov character [28-30 . Within the zero-range 
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the treatment by Endo et al. does not seem to yield an 
analogous physical picture. Additionally we speculate 
that while Ref. [T^] employs two different models for the 
momentum cutoff A c , it docs not explore the entirety of 
the (Ro,S(Rq)) parameter space. In the future it will be 
interesting to investigate the precise connection between 
the formulations in the coordinate and momentum spaces 
by, e.g., comparing the wavef unctions. Such a compari- 
son is needed to check if the above correspondencies are 
correct. 



CONCLUSION 



FIG. 7. (Color online) Circles, diamonds and squares show 
the off-diagonal coupling ft/ E2b used in the two-state model, 
Eq. {Uj ), as a function of k for R /a s = 0.0001,0.0003 and 
0.001, respectively, on a logarithmic scale. We used E u — E u ,j 
in the two-state model. The dependence of /?/ E2b on the value 
of Ro increases with increasing k. 



framework employed here, the exact number and energy 
of the Efimov trimers supported depends on the short- 
range hyperradial boundary condition. If the hyperradial 
boundary condition is fixed and k is varied, the energy 
spectrum changes smoothly from (i) deviating from the 
universal spectrum only in a small region around S c to (ii) 
deviating from the universal spectrum for a fairly large 
range of 6 to (hi) supporting three-body states with Efi- 
mov character for k > 13.606. Thus, the deviations from 
universality discussed in this paper can be interpreted as 
connecting the universal states predicted by Kartavtsev 
and Malykh and the Efimov trimers. 

Our findings are in qualitative agreement with the re- 
sults reported in Ref. [T2]. In that work, the authors 
imposed a momentum cutoff A c in the Skorniakov-Ter- 
Martirosian equation and used the quantity £?~ E ^ b 
(i = 1 or 2) to determine the "boundary" between the 
universal trimers predicted by Kartavtsev and Malykh 
and the crossover trimers in the (k, (A c as)^ 1 ) parame- 
ter space. The momentum cutoff was introduced in two 
ways, using a sharp and a Gaussian cutoff. We speculate 
that in our formulation a change in the three-body phase 
and/or hyperradius Rq corresponds to a change in A c . A 
crossover trimer, in turn, corresponds to a trimer whose 
energy deviates appreciably from E^^. While the devi- 
ations from E u ^ are, in our formulation, linked to the 
coupling of the universal state to a non-universal state, 



In this paper we studied a system of two identical 
heavy fermions of mass M and a light particle of mass m 
with zero-range two-body interspecies interactions. In 
particular, we looked at deviations of the three-body 
bound state energies from the universal energies E u ^i for 
k > Ki and E Ut i and E u> 2 for K2 < k < 13.606 as a func- 
tion of the hyperradial short-range boundary condition. 
We imposed a short-range phase 5(Rq) using a logarith- 
mic derivative boundary condition at various hyperradii 
Rq. This parameterization allowed us to explore the full 
range of possible short-range boundary conditions. 

We found that (i) for 8(Rq) — ir/2 the universal states 
with energies E u i and E u 2 , predicted by Kartavtsev and 
Malykh [TT], are recovered; (ii) the three-body states de- 
viate from universality in the vicinity of a three-body res- 
onance located at the short-range phase S c (Rq), at which 
the non-universal state is first bound; (hi) the deviations 
from universality increase with increasing mass ratio k 
(at fixed Rq) and with increasing Rq (at fixed k); (iv) 
QDT accurately predicts the values of S c (Rq); (v) a two- 
state model quantitatively describes the behavior of the 
universal and non-universal states as a function of <5(i?o); 
and, finally, (vi) the non-universal bound state exists for 
k < Ki even though universal bound states are not sup- 
ported in this regime. 
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